Entanglement of Atomic Qubits using an Optical Frequency Comb 
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We demonstrate the use of an optical frequency comb to coherently control and entangle atomic 
qubits. A train of off-resonant ultrafast laser pulses is used to efficiently and coherently transfer 
population between electronic and vibrational states of trapped atomic ions and implement an 
entangling quantum logic gate with high fidelity. This technique can be extended to the high field 
regime where operations can be performed faster than the trap frequency. This general approach 
can be applied to more complex quantum systems, such as large collections of interacting atoms or 
molecules. 
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The optical frequency comb generated from an ultra- 
fast laser pulse train has revolutionized optical frequency 
metrology [T]-[4] and is now playing an important role in 
high resolution spectroscopy [5] . The spectral purity yet 
large bandwidth of optical frequency combs also provides 
a means for the precise control of generic quantum sys- 
tems, with examples such as the quantum control of mul- 
tilevel atomic systems [6l [7] , laser cooling of molecules or 
exotic atomic species O [9] , and quantum state engineer- 
ing of particular rovibronic states in molecules [lOl [11] . 
The optical frequency comb may become a crucial com- 
ponent in the field of quantum information science, where 
complex multilevel quantum systems must be controlled 
with great precision [T2^. 

In this Letter, we report the use of an optical frequency 
comb generated from an ultrafast mode-locked laser to 
efficiently control and faithfully entangle two trapped 
atomic ion qubits. The optical pulse train drives stimu- 
lated Raman transitions between hyperfine levels p!3lfT4]. 
accompanied by qubit state dependent momentum kicks 
|15 . The coherent accumulation of these pulses generates 
particular quantum gate operations that are controlled 
through the phase relationship between successive pulses. 
This precise spectral control of the process along with the 
large optical bandwidth required for bridging the qubit 
frequency splitting forms a simple method for control- 
ling both the internal electronic and external motional 
states of trapped ion qubits, and may be extended to 
most atomic species. This same approach can be applied 
to control larger trapped ion crystals with more advanced 
pulse-shaping techniques, and can also be extended to a 
strong pulse regime where only a few high power pulses 
are needed for fast quantum gate operations in trapped 
ions [ISHni. 

High fidelity qubit operations through Raman transi- 
tions are typically achieved by phase-locking frequency 
components separated by the energy difference of the 
qubit states. This is traditionally accomplished in a 
bottom-up type of approach where either two monochro- 



matic lasers are phase-locked or a single cw laser is modu- 
lated by an acousto-optic (AO) or an electro-optic (EO) 
modulator. However, the technical demands of phase- 
locked lasers and the limited bandwidths of the modu- 
lators hinder their application to experiments. Here we 
exploit the large bandwidth of ultrafast laser pulses in 
a simple top-down approach toward bridging large fre- 
quency gaps and controlling complex atomic systems. By 
starting with the broad bandwidth of an ultrafast laser 
pulse, a spectral landscape can be sculpted by interfer- 
ence from sequential pulses, pulse shaping and frequency 
shifting. In this paper, we start with a picosecond pulse 
and, through the application of many pulses, generate a 
frequency comb that drives Raman transitions by stim- 
ulating absorption from one comb tooth and stimulating 
emission into another comb tooth as depicted in Fig. [l] 
Because this process only relies on the frequency differ- 
ence between comb teeth, their absolute position is ir- 
relevant and the carrier-envelope phase does not need to 
be locked [2] . As an example of how this new technique 
promises to ease experimental complexities, the control 
of metastable-state qubits separated by a terahertz was 
recently achieved using cw lasers that are phase-locked 
through a frequency comb [18 , but might be controlled 
directly with a 100 fs Ti:sapph pulsed laser. 

At a fixed point in space, an idealized train of laser 
pulses has a time-dependent electric field that can be 
written as 

N 

E{t) = Y,f{t-nT)e'"^\ (1) 

n=l 

where f(t) is the pulse envelope, T is the time between 
successive pulses (repetition rate z^r = 1/T), N is the 
number of pulses in the train and ujc is the carrier fre- 
quency of the pulse. For simplicity, any pulse-to-pulse 
optical phase shift is ignored since the offset frequency 
in the comb is unimportant. The Fourier transform 
of Eq. ([T]) defines a frequency comb characterized by 
an envelope fioo) = ^[f{t)] centered around the opti- 
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FIG. 1: The Stokes Raman process driven by frequency combs 
is shown here schematicahy. An atom starting in the |^) state 
can be excited to a virtual level by absorbing a photon from 
the blue comb and then driven to the |t) state by emitting a 
photon into the red comb. Although drawn here as two dif- 
ferent combs, if the pulsed laser's repetition rate or one of its 
harmonics is in resonance with the hyperfine frequency, the 
absorption and emission can both be stimulated by the same 
frequency comb. Because of the even spacing of the frequency 
comb, all of the comb teeth contribute through different vir- 
tual states which result in indistinguishable paths and add 
constructively. 



cal frequency uOc and teeth separated by whose in- 
dividual widths scale like ~ vr/N. The Raman res- 
onance condition will be satisfied when a harmonic of 
the repetition rate is equal to the hyperfine qubit split- 
ting cjo, implying that the parameter q = ujq/27:vr^ 
is an integer. To demonstrate coherent control with a 
pulse train, ^^^Yb^ ions confined in a linear Paul trap 
are used to encode qubits in the ^6'i/2 hyperfine clock 
states \F = 0,mF = 0) = ||) and \F = l,mF = 0) = |t), 
having hyperfine splitting uoq /2ti — 12.6428 GHz. For 
state preparation and detection we use standard Doppler 
cooling, optical pumping, and state-dependent fluores- 
cence methods on the 811 THz ^5'i/2 ^ ^^1/2 electronic 
transition The frequency comb is produced by a 

frequency-doubled mode- locked Ti: Sapphire laser at a 
carrier frequency of 802 THz, detuned by A/27r = 9 THz 
from the electronic transition. The repetition rate of the 
laser is vr = 80.78 MHz, with each pulse having a du- 
ration of r ^ 1 psec. The repetition rate is phase-locked 
to a stable microwave oscillator as shown in Fig. [2j pro- 
viding a ratio of hyperfine splitting to comb spacing of 
q = 156.5. An EG pulse picker is used to allow the pas- 
sage of one out of every n pulses, decreasing the comb 
spacing by a factor of n and permitting integral values 
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FIG. 2: Schematic of the experimental setup showing 
the paths of the pulse trains emitted by a mode-locked 
Ti:Sapphire (Ti:Sapph) laser, where the optical pulses are fre- 
quency shifted by AOs. Single qubit rotations only require a 
single pulse train, but to address the motional modes the pulse 
train is split into two and sent through AOs to tune the rela- 
tive offset of the two combs. We lock the repetition rate (i^r) 
by first detecting ur with a photodetector (PD). The output 
of the PD is an RF frequency comb spaced by ur. We band- 
pass filter (BP) the RF comb at 12.685 GHz and then mix the 
signal with a local oscillator (LO). The output of the mixer is 
sent into a feedback loop (PID) which stabilizes vr by means 
of a piezo mounted on one of the laser cavity mirrors. When 
locked, ur is stable to within 1 Hz for more than an hour. 
As an alternative, instead of locking the repetition rate of the 
pulsed laser, an error signal could be sent to one of the AOs 
to use the relative offset of the two combs to compensate for 
a change in the comb spacing. 



of q. As shown in Fig. [sj when n = 2 {q = 313 and 
vr = 40.39 MHz), application of the pulse train drives 
oscillations between the qubit states of a single ion. How- 
ever, when n = 3 {q = 469.5 and ur = 26.93 MHz), the 
qubit does not evolve. 

The Rabi frequency of these oscillations can be esti- 
mated by considering the Hamiltonian resulting from an 
infinite train of pulses. After adiabatically eliminating 
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FIG. 3: After Doppler cooling and optical pumping to the \l) 
state, a single pulse train is directed onto the ion. When the 
ratio of qubit splitting to pulse repetition rate, q, is an integer, 
pairs of comb teeth can drive Raman transitions as shown by 
the blue circular data points. However, if the q parameter is 
a half integer, the qubit remains in the initial state as shown 
by the red square data points. 



3 



the excited state and performing the rotating- wave 

approximation, the resonant Rabi frequency of Raman 
transitions between the qubit states is given by a sum 
over ah spectral components of the comb teeth as indi- 
cated in Fig. [l](n = 1): 



UJqT 



(2) 



where /i is the dipole matrix element between the ground 
and excited electronic states, E]^ = i^Rfi^Trkun)^ and 
q is an integer. In the approximate expression above, 
the sum is replaced by an integral and each pulse is de- 
scribed by f{t) = ^/7T/2Eosech{7^t/r) with r <^T^ where 
^0 = {^R^)\fJ^Eo\'^ /A = 57^/2 A is the time- averaged res- 
onant Rabi frequency of the pulse train and s = I/Isat is 
the average intensity / = iyjiceo/2 J dt\f{t)\'^ scaled to the 
^ saturation intensity. Note the net transi- 

tion rate is suppressed unless the single-pulse bandwidth 
is large compared to the hyperfine frequency {ujqt <C 1), 
in which case Q ^ Qq. In our experiments, c^qt ~ 0.08. 



For Isat = 0.15 W/cm^, the data shown in Fig. [sjis con- 
sistent with an average intensity I ~ 500 W/cnr. 



In order to entangle multiple ions, we first address the 
motion of the ion by resolving motional sideband tran- 
sitions. As depicted in Fig. |2| the pulse train is split 
into two perpendicular beams with wavevector difference 
k along the x— direction of motion. Their polarizations 
are mutually orthogonal to each other and to a weak 
magnetic field that defines the quantization axis [20] . 
We control the spectral beatnotes between the combs by 
sending both beams through AO modulators (driven at 
frequencies ui and 1^2)^ imparting a net offset frequency 
of A(jj/27T = vi — V2 between the combs. For instance, 
in order to drive the first upper/lower sideband transi- 
tion we set \2T:jvji + Ac^l = uoq ^ uOf, with j an integer 
and uot the trap frequency. In order to see how the side- 
bands are spectrally resolved, we consider the following 
Hamiltonian of a single ion and single mode of harmonic 
motion interacting with the Raman pulse train: 



J 



Y 



(3) 



where 6p = is the change in the Bloch angle due to 
a single pulse, cf^ is the Pauli-z operator, cf± are raising 
and lowering operators, x is the x— position operator of 
the trapped ion, o) and a are the raising and lowering 
operators of the x— mode of harmonic motion and the q 
parameter has been assumed to not be an integer or half- 
integer. In the interaction picture, the evolution operator 
after N pulses is given by ^ where 



V = exp [-iHoT] exp 



-iO, 
~2 



(4) 



and Hq = ujto^a + 1/2 (cjq + Acjja^. The time evolution 
operator is given by, 

N-l 

= exp[-zi7oiVT](/ - E + ^(^p)) (^) 

n=0 

i(u;o+Au;)nT -^^iujtnT^ 



L>(z7?e*^*^' ) + h.c. , 



(6) 



where D{a) = exp[aaT — a* a] is the harmonic oscil- 
lator displacement operator in phase space, and r] = 
k^h/2mujt is the Lamb-Dicke parameter. In the 
Lamb-Dicke regime, r]^ (a^a) + 1 <C 1, we can write 
Diirje'^'"^^) ^ 1 + z7/(e^^*^^at + e-*^*^^a) turning the 
sum in Eq. ([5| into a geometric series. If, for exam- 
ple, the offset frequency between the combs Auj is tuned 



to satisfy the resonance condition for the red sideband, 
i9r = {ojQ + — ^t)T = 27rj, where j is an integer, then 
the sum in Eq. ([5| is approximately given by. 



N-l 
n=0 



.^siniVV2gitf,(^- 
sin'i^^/2 



•i)/2 



cr+a + h.c. (7) 



The coefficient in Eq. ^ is the same as the field ampli- 
tude created by a diffraction grating of N slits, whose 
narrow peaks have an amplitude equal to N. In the 
limit uJtT <C 1, the other terms in Eq. ([5| that drive the 
carrier and other sideband transitions can be neglected 
when N {ujtTr])~^ . This is analogous to the de- 
structive interference of amplitudes away from the bright 
peaks in a diffraction grating. For Ut/27T = 1.64 MHz, 
T = 12.4 ns and r] = 0.1, the sidebands are well-resolved 
when TV > 80. 

For many applications in quantum information, the 
motional modes of the ion must be cooled and initial- 
ized to a nearly pure state. Fig. [4] shows that the pulsed 
laser can also be used to carry out the standard tech- 
niques of sideband cooling [20] to prepare the ion in the 
motional ground state with near unit fidelity. The set-up 
also easily lends itself to implementing a two-qubit entan- 
gling gate by applying two fields whose frequencies are 
symmetrically detuned from the red and blue sidebands 
[2T| [22] . By simultaneously applying two modulation fre- 
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FIG. 4: (a) Using a Raman probe duration of 80/is, (N ^ 
6500), a frequency scan of AOl shows the resolved carrier and 
motional sideband transitions of a single trapped ion. The 
transitions are labeled, {Aux, Any), to indicate the change in 
the number of phonons in the two transverse modes that ac- 
company a spin flip. The x and y mode splitting is controlled 
by applying biasing voltages to the trap electrodes. Unla- 
beled peaks show higher order sideband transitions and tran- 
sitions to other Zeeman levels due to imperfect polarization of 
the Raman beams, (b) Ground state cooling of the motional 
modes via a train of phase-coherent ultra-fast pulses. The 
red open-circle data points show that after Doppler cooling 
and optical pumping, both the red and blue sidebands are 
easily driven. The blue filled-circle data points show that af- 
ter sideband cooling, the ion is close to the motional ground 
state, {nx,y < 0.03), as evidenced by the suppression of the 
red-sideband transition. 



quencies to one of the comb AO frequency shifters, we 
create two combs in one of the beams. When these combs 
are tuned to drive the red and blue sidebands (in conjunc- 
tion with the third frequency comb in the other beam), 
the ion experiences a spin-dependent force in a rotated 
basis as described in Ref. [23]. Ideally, when the fields 
are detuned from the sidebands by an equal and opposite 
amount 6 = 2/^1], a decoupling of the motion and spin oc- 
curs at gate time tg = 27r/(5, and the spin state evolve to 
the maximally entangled state |x) = + e*"^ In 
the experiment, tg = 108/is (TV ~ 8700 pulses). 

The entanglement is verified by the measurement of 
a fidelity-based entanglement witness [24 . The fidelity 
of the state p with respect to |x) is found by measur- 
ing the populations p^^,^^ and Pn^n and scanning the 
measurement phase angle (j) in the parity signal 11(0) = 
Tr[ai^^ai^^R{(t)f^pR^cl)f^], where i?(0) is a 7r/2 rota- 




FIG. 5: The parity oscillation that is used to calculate the 
fidelity of the spin state of two ions with respect to the max- 
imally entangled state |x) after performing the entangling 
gate. The phase (j) of the analyzing pulse is scanned by chang- 
ing the relative phase of the rotation pulses. The offset and 
lack of full contrast in the parity signal can be attributed to 
state detection errors. 



tion on the Bloch sphere with phase [25] . The contrast 
of the parity signal, lie, is used to calculate the fidelity 
J- = {p^^^^^ -f- pii^^i)/2 + nc/4. The measured popula- 
tions of III) and Itt) together with the data shown in 
Fig. [5] yield a fidelity T = 0.86 ± 0.03, thereby signaling 
that the two ions are entangled after the application of 
the pulse train. 

We have demonstrated full control and entanglement of 
two atomic qubits using an optical frequency comb. This 
work represents a significant simplification over current 
methods for optical control of trapped ion qubits, and 
also points the way toward future advances with higher 
power laser pulses. For example, such pulses allow much 
larger detunings from resonance and a suppression of de- 
coherence from spontaneous emission [26 while main- 
taining gate speed. When only a few high-power pulses 
are used in ways similar to the experiment, it also be- 
comes possible to suppress sources of motional decoher- 
ence through the use of fast entanglement schemes [16 . 
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